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Based on an effective model of a doped antiferromagnetic Mott insulator, we show that a doped 
hole will induce a dipole-like spin configuration in a spin ordered phase at low doping. The kinetic 
energy of doped holes is severely frustrated and a hole-dipole object is actually localized or self- 
trapped in space. Without a balance from the kinetic energy, the long-range dipole-dipole interaction 
between doped holes will dominate the low-energy physics, leading to an inhomogeneity instability 
as hole-dipoles collapse into stripes. Both antiphase metallic stripes of quarter-filling and antiphase 
insulating stripes along a diagonal direction are discussed as composed of hole-dipoles as elementary 
building blocks. Stripe melting and competing phases are also discussed. 
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I. INTRODUCTION 

The theoretical challenge of doped antiferromagnetic 
(AF) Mott insulator systems comes from their nature 
of strong correlations. It makes the relevant models like 
the Hubbard and t — J Hamiltonians highly difficult to 
deal with. In such systems, the kinetic enerry of charge 
carriers usually gets strongly frustrated 2, 3] and there 
is no single dominant process based on which a pcrtur- 
bative approach can be constructed. This opens door for 
various competing orders at low doping and makes the 
phase diagram complex 4]. 

So far there still lacks a single microscopic description 
which can unify various competing phenomena known 
from both experiments in cuprates and theoretical con- 
siderations based on doped antiferromagnets. Given the 
fact that exact solutions of two-dimensional (2D) Hub- 
bard and t — J models are absent, a simplified and 
tractable model of doped AF Mott insulators is thus 
highly desirable in order to have a framework for a sys- 
tematic study. 

Such a type of models should incorporate three most 
essential ingredients of a doped AF Mott insulator. 
Namely, it should properly describes AF correlations at 
half-filling, including both long-range (low-energy) and 
short-range (high-energy) ones; Its Hilbert space should 
be restricted due to the presence of an upper Hubbard 
band at finite doping; Finally, it should explicitly embody 
essential frustration effects introduced by hole hopping in 
an antiferromagnet. Furthermore, compared to the orig- 
inal Hubbard and t — J models, it ought to be simplified 
enough to allow a tractable treatment. 

In this paper, we study an effective model @ of 
a doped AF-Mott insulator which satisfies the above 
criteria. This model well describes AF correlations at 
half-filling and precisely incorporates the so-called phase 
string effect at a finite doping. The latter has been identi- 
fied jj as a crucial dynamic frustration effect introduced 
by holes in the t — J Hamiltonian. Here a phase string 
depicts a string of displaced transverse spins along a hole- 



hopping path, characterized by a sequence of signs: 



(+1) X (-1) X (-1) X ... (1) 



as illustrated by Fig. 1(a). It can be rigorously shown 
Q that such transverse spin mismatches cannot be self- 
healed simultaneously with the displaced spins in the lon- 
gitudinal (quantization) direction [see Fig. 1(b)], through 
the superexchange term of the t — J Hamiltonian. Con- 
sequently, after a hole moves through a closed path back 
to its origin and the displaced spins by hole hopping 
are restored in the quantization direction (i-axis), the 
wavefunction should always pick up an additional Berry- 
phase-like sign given by which is highly nontrivial in 
a general case. 

Such a working model may be called the phase string 
model (section II). Based on it, we will show that in a spin 
ordered phase at low doping, a doped hole will generally 
induce a dipole-like spin configuration (section HI). The 
kinetic energy of doped holes is severely frustrated by the 
phase string effect such that a hole-dipole is localized or 
self-trapped in space with an infinite mass. Without the 
balance from the kinetic energy, the long-range dipole- 
dipole interaction between doped holes will dominate the 
low-energy physics. It will lead to an inhomogeneous in- 
stability as hole-dipoles collapse into stripes (section IV). 
These anti-phase stripes can be understood as that holes 
are threaded together by the phase strings. Both metal- 
lic stripes of anti-phase and quarter-filling and insulating 
stripe along the diagonal direction are discussed based 
on the elementary hole-dipoles. The stripe phase at low 
doping is highly competitive to homogeneous phases pre- 
viously studied 01 in the same model, which include the 
AF long-range ordered state, spin glass state, and super- 
conducting state. Therefore, a rich array of phases are 
revealed in the phase string model, where AF correlations 
and phase string effect serve as basic driving forces. 
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FIG. 1; Spin mismatches created by a hole hopping in a Neel 
order: (a) in transverse directions (phase string) and (b) in 
longitudinal direction. Two dotted circles in (a) indicate the 
sigularity of the starting and ending points in a holon path 
(see text). 



II. AN EFFECTIVE MODEL OF A DOPED 
ANTIFERROMAGNET 

The t — J Hamiltonian is one of the simplest nontrivial 
models of doped Mott insulator/doped antiferromagnets 
0. Based on the exact phase string formulation j^] of 
the t — J Hamiltonian, an effective low-energy model can 
be derived, known as the phase string model 0. 

The presence of short-range AF correlations provides a 



natural "ultraviolet" cutoff in deducing such an effective 
model from the t — J Hamiltonian. So the local AF cor- 
relations, characterized by the bosonic RVB pairing 
underpin this effective theory. It differs from the slave- 
boson mean-field theories by emphasizing the under- 
lying AF correlations; It differs from the slave- fermion 
mean-field theories [Tfl llll | by precisely dealing with a 
singular hopping effect in the t — J model, i.e., the phase 
string effect mentioned in Introduction. 

Such an effective model itself defines a class of doped 
Mott-insulator/doped antiferromagnet. It has been 
shown to possess the following properties: (I) It well 
describes antiferromagnetism at half-filling, including 
both long-wavelength, low-energy and short-range, high- 
energy (high-temperature) physics; (III) It has a d-wave 
superconducting ground state beyond a finite doping con- 
centration; (II) It characterizes the correct restricted 
Hilbert space for charge and spin degrees of freedom of a 
doped Mott insulator, when holes are introduced. 

The effective Hamiltonian can be written [H, |^ as 
Hst 



ring 



Hh + Hg , where 



Hh = -th^(e'^-^-<)/iI/i,+F.c. (2) 
= -J,^ (e-^?.)6L6]_, + i7.c. (3) 

with hi and bia- denoting bosonic holon and spinon oper- 
ators, respectively. When the holon number is zero, Hs 
is equivalent to the Schwinger-boson mean-field Hamil- 
tonian J which well describes AF correlations at half- 
filling. The AF long range order is realized by the spinon 
Bose condensation, i.e., (bi^) ^ ai zero temperature. 
At larger doping, when holons are Bose condensed, the 
ground state becomes a (d-wave) superconductor 

A unique feature of this effective model is the presence 
of link fields in Eqs. ^ and They are not free gauge 
fields as they satisfy the topological constraints 



E4 = ±-E-'' 



and 



E 



A 



(4) 



(5) 



for a closed loop c, in which and are holon and 
spinon number operators, respectively. In addition to 
Al^ , (jJ^l^ in Eq. Q describes a uniform tt flux per plaque- 

tte (X^n 4>% = i"")- This model has a U(l)xU(l) gauge 
symmetry as Hh and Hg are invariant under the gauge 
transformations: 



hj exp{iipj), Alj Alj -t- i{ipi - ipj), (6) 



and 



6j^exp(icr0j)> ^ii 



(7) 
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Here the holons and spinons are basic building blocks 
that properly characterize the restricted Hilbert space of 
a doped Mott insulator. But it is important to note that 
they, as independent entities, are not necessarily the true 
low-lying elementary excitations of the system. In this 
model, the electron operator is made of holon and spinon 
operators in the following form (5j 



= h]b,,^ 



(8) 



where Qia- is a topological phase, which ensures the 
fermionic statistics of Cio-, as defined by 



with 



$^ = ^Imln(z, 



and 



^0 E- 



Iniln(zi — zi)n^ 



(9) 



(10) 



(11) 



Thus, going from the holon-spinon representation to the 
electron representation, topological (vortexlike) phases, 
$^ and are involved, which implies that the holon 
and spinon objects defined in the phase string model are 
topological entities themselves, as will soon become clear 
below. 



III. DOPED HOLES: DIPOLES IN SPIN 
ORDERED PHASE 

How doped holes propagate in an antiferromagnet and 
how background spins are influenced by the hole motion 
are two central issues in the study of doped antiferro- 
magnets. In the present section, we will focus on these 
issues in a spin ordered phase at small doping, based on 
the above phase string model. 



Holons as merons 



The spin-flip operator can be expressed by 



(12) 



according to Eq. |(SJ). The presence of the vortex phase 
$^ in Eq. H12|l clearly shows that a holon is a topological 
object that affects spins nonlocally. 

At half-fllling where $^ = 0, the Bose condensation < 
bia of spinons gives rise to an AF long-range order 
of spins lying in the x-y plane: (5*^^) = (^1)' < b]^ >< 
bii >. Through each holon will induce a "meron" 
spin twist, which holds true as long as < 6^0- >^ 0. Such 
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FIG. 2: A meron-like spin configuration around a holon. 



a spin vortex configuration (for a holon at the origin) 
may be characterized by the unit vector 



as follows: 



n, (x(™l)* (S,), 



n,= Vil I 



(13) 



with r| = -t- y1. It is schematically shown in Fig. 2. 

To evaluate the potential energy of a holon (meron), 
we consider the spinon condensed part in [Eq. (O] : 

AS, ~ -JsPIo^ [cos A% - 1] + c.c. (14) 

in which p^a^ =< 6^ >^(a is the lattice constant) and 
a reference energy corresponding to A^- = has been 

subtracted. Explicitly J^, reads 



4 = i^Imln 



Z^ - h 
Zj - Zl 



(15) 



and then one has 



(16) 



where i?" 
energy 



J2i n'lS™' with a single holon-meron self- 



JsPlaHn{L/a), 



(17) 
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where L is the size of the sample. Namely, a single holon 
as a meron in the spinon condensed phase will generally 
cost a logarithmically divergent energy. 

Furthermore, the meron-meron interaction is given by 



Kl' 



and 



V™' 
II' 



JsPla 



2 S 



Im In ■ 



(18) 



— Im In 



Zl 



-nJspla^ In 



(19) 



Here two holons will always repulse each other loga- 
rithmically, as they are merons of the same topological 
charge. 

We note that there have been proposals of hole-merons 
in literature ^13^ ,1_4] . They usually involve both merons 
and anti-merons . Two holes may thus be bound like a 
meron-anti-meron (vortex-antivortex) pair. By contrast, 
in the present case holons are always repulsive to each 
other as they carry the same topological charge. The 
underlying physical reason is that this is in favor of the 
kinetic energy for a finite concentration of holons, which 
become identical bosons after explicitly incorporating the 
phase string effect in the slave-fermion formalism. Conse- 
quently, superconductivity in the present doped AF-Mott 
insulator originates from the Bose condensation of 
holons, instead of the pairing of them. Here one should 
distinguish the holon pairing from the pairing order pa- 
rameter of electrons. The latter still characterizes the su- 
perconducting state in the present case where spins are 
in RVB pairing, even though the holon pairing is absent. 

Intuitively one might envisage that a phase string in- 
duced by one hole be "erased" by an another hole accom- 
panying it. This may look like in favor of hole pairing. 
But such a picture is actually more in favor of stripes 
in which many holes precisely follow each other's paths 
to eliminate the phase string effect. The basic reason is 
that the phase string effect of ^ is not simply a "geo- 
metric" effect: it also involves spin dynamics and is very 
singular. Thus, unless one hole tightly follows the other 
along the same path, like in the stripe case, the singu- 
lar dynamic phase string effect cannot be completely re- 
moved by two holes simply staying close to each other. 
In fact, previously Trugman 0,3 has already found the 
presence of frustrations for a hole pair's motion in an an- 
tiferromagnet, due to a quantum effect originating from 
the fermionic character of the background spins. This 
effect is actually included in the phase string effect in 
the present approach. Finally, we note that technically, 
a holon is just a building block in the phase string model 
and the possibility that some holons may be anti-merons 
can be, in principle, constructed as composite objects in 
the present framework. But such a construction is usu- 
ally not energetically favorable at finite doping and will 
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FIG. 3: A hole-dipole is a confined object composed of a 
holon-meron and an anti-meron at two poles: B and A, con- 
nected by a branch-cut with the spatial separation e as a 
dipole moment. 



be discussed elsewhere. The pairing of two doped holes 
has been discussed fll| in the t — Jz model based on the 
"string" picture before. 



B. Holon confinement: Doped holes as dipoles 

To remove such an energy divergence [Ea. (|17|l ] asso- 
ciated with a meron spin configuration, each holon will 
have to "nucleate" an anti-meron from the vacuum. Be- 
ing bound to the latter, a dipole composite with a finite 
energy will be formed. This elementary object will carry 
a charge +e and replace the holon to represent a true 
stable hole in the spin ordered phase. By contrast, free 
holons disappear from the finite-energy spectrum. We 
can also call this phase as the holon confined phase. 



1. Formation of dipoles 



An anti-meron configuration is defined as a spin twist 



bia exp 



(20) 



where 'd^ = Imln(zi — z°). Here denotes the position 
(center) of an anti-meron labeled by k. Correspondingly 



{S+) ^ (-1)' (bl^) (fe.x) exp [z^': - iQ,] (21) 

where Hi = (summing over all induced anti- 

merons). If z^'s are close to the positions of correspond- 
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ing holons, then cancellation occurs in $^ — ^i- Conse- 
quently, the logarithmic divergence is removed from Hs 
(see below). 

Define nf + inf = e**^""*""^*, with the unit vector no = 
(cos (jjo, sin(/)o) as the magnetization direction at infinity. 
For a bound pair of meron (holon)-anti-meron centered at 
the origin, labeled by the index k, one obtains according 
to (EH 



= Im In 



Zfc/2 



(22) 



with Zk = + ie\ . Here denotes the spatial displace- 
ment of the meron and anti-meron centered at re- 
spectively. At |ri| >> \ek\, 

,k _ (z X efc) 

- — nn? — • (23) 











•—J 





Type A dipole 

(a) 



Type B dipole 

(b) 



Thus, each pair of holons and anti-merons forms a com- 
posite, as shown in Fig. 3, that has a dipolar spin con- 
figuration at a sufficiently large distance. 

In Hs, one may easily check that uJij = (il^ — J7j)/2 can- 
cels out A'lj and thus removes the aforementioned energy 
divergence. Here the energy cost Eq. (|14f) is replaced by 

AEs = -2Jspla^ [co^ i^^3 - '^y) ~ 1] 
J,p^a2^(4-c.,,)' 

k<k' 

in which the single dipole energy 



cd _ 
^k — 



and dipole-dipole interaction energy 



d-d 
kk' 



In the continuum limit, 
according to Eq. (j^ . 



^ (r. 



(25) 

• (26) 
V0''" and, 



Vet)" 



zx(r-efe/2) zx(r + efc/2) 



|r-efc/2| 



-e,/2r 



(27) 



Then the self-energy of a dipole can be evaluated by Sjj: 



J (fir(Vcf)'^) ■ In contrast to the logarithmically 



two poles at ±^ of the hole dipole, and £^ becomes finite 



divergent meron-energy, there is a branch-cut connecting 
■ oft] 

4 

~ 7rJ,p!aMn ' (28) 
a 

at lefcl > a. 



^fc| 

r-efc/2|'|r + efc/2|2 



FIG. 4: Two typical hole-dipoles of minimal size: (a) Type A 
dipole with a holon and a spin sitting at two poles at nearest- 
neighboring lattice sites; (b) Type B dipole with one pole at 
the center of a plaquette. 



2. Type A and Type B dipoles 

Let us first point out that generally there are two pos- 
sible anti-meron configurations as defined by Eq. I|20|l . 
Type A: the core sits on a lattice site /, i.e., 2;° = z;. In 
this case, there will be always a spin 1/2 sitting at the 
core of the anti-meron [Fig. 3], according to the single- 
occupancy constraint. Type B: z^l is located inside a 
lattice plaquette, instead of a lattice site. Figs. 4(a) and 
(b) illustrates both cases at minimal sizes. 

Type A dipole. Besides the energy consideration given 
above, it is important to make sure the self-consistency 
of the theory, that is, the anti-meron configuration Eq. 
(Pn|l should not violate the single-valueness of the elec- 
tron operator (the spin operator in Eq. H21|l is obviously 
single- valued). 

Consider that a winding of the coordinate Zi around 
2:° which is located at a lattice site. A phase ±7r (cor- 
responding to cr = ±1) will be acquired on the r.h.s. of 
Eq. (pn|l . In order to ensure the single-valueness of Cia, 
the holon field has to be transformed in response to Eq. 
(Pn|l . We may introduce 



In (2:.-z0) 



(29) 



with the spin index ak = ±1 denoting the spin at the 
site z^. Then the total change in the electron operator 

be 



(for a winding of the coordinate Zi around z^) will 

„i((Tir-afc7r) _ ^ g-j-jj^j satisfying the single-valueness 



condition. 



6 



The physical imphcation of H29|l is as foUows. It will 
change in Eq. |(2Jl by 

A^-.4^A^-^Imln|i^. (30) 

It simply removes the contribution of the spinon, 
(jn\^ = ak, at the site of zi = z'^ from Afj, i.e., 

5]i^=^A^-afc^, (31) 

c c 

where c is a loop encircling z{?. By incorporating the phase 
on the r.h.s. of Eq. (^i^ into 6^^ of Eq. (jHJ, one may 
also see that the contribution from the spinon at the site 
of z1 is eliminated from <i>^ and thus from Cicr- It means 
that the spin at z^ will no longer exert the nonlocal effect, 
embedded in a spinon, through Af^ in and Qia- in Cia- ■ 
Hence, to maintain the single-valueness condition after 
the creation of an anti-meron, the spin 1/2 at its center 
site z^ will lose its spinon character and looks like a local 
"defect" spin. 

Such a "defect" spin, sitting at the core of an anti- 
meron, will be decoupled from the environmental spins 
and simply behave like a free 5 = 1/2 spin. This can 
be straightforwardly checked. By noting that the trans- 
formation Eq. (pn|) excludes the case Zi = z° due to the 
uncertainty of i?*^, the phases of those terms in Hg that 
involve the links attached to the site 2°, will become un- 
certain after the transformation. Specifically, these terms 
connected to z^ will all vanish when one averages over i9o 
after the replacement 'd'^ 'd'^ + i?o in Eq. ipH]! for 
Zi 7^ z^, while the rest of Hg remains invariant. Thus, 
the spin 1/2 at z'^ is indeed decoupled from the environ- 
mental spins. 

Type B dipole. Now consider the case that z° is located 
at a center of a plaquette [see Fig. 4(b)]. There is no 
longer a lattice spin at and the holon transformation 
Eq. (|29|l may be replaced by the form 

h^^hle^'^'l (32) 

One can check that Cia- still remains single-valued for a 
winding of the coordinate Zi around z'f^,. Based on Eq. 

, the holon Hamiltonian Hh is changed by Afj Af^ , 
with 

if^.^A^iilmln^i^. (33) 

Physically it means that an additional ±7r flux penetrates 
through the plaquette where resides: 

c c 

with c as a loop encircling the given plaquette. 

In contrast to Type A, it is not obvious that there is a 
free spin 1/2 trapped around the plaquette of the core of 



the anti-meron. In Type B, the cost of the superexchange 
energy is mainly due to the spin twist around the plaque- 
tte zi which is less dramatic than in Type A, where the 
presence of an isolated spin 1/2 at z° means the loss of the 
superexchange energy for four nearest-neighboring links 
connected to it. But in Type B, there are additional tt 
flux quanta [see Eq. H34|l ] threading through the plaque- 
ttes, where the cores of anti-merons are located, seen by 
holons. So the hopping energy of holons in Type B is dif- 
ferent from that of a Type A dipole. Generally, Type A 
dipoles are expected to be more stable at higher doping, 
while Type B dipoles should be energetically favorable at 
low doping. 



3. Localization of hole dipoles and phase string effect 

The fact that a doped hole will induce a dipolar spin 
configuration was first pointed out by Shraiman and Sig- 
gia |la| based on a semi-classical treatment of the t ~ J 
model. In Shraiman-Siggia's picture, the hole is sitting at 
the center of the dipole and such a hole-dipole is mobile in 
space. This mobility may be easily understood if one con- 
siders the translational invariance and the fact that the 
hopping overlap integral for two hole-dipolcs with their 
centers displaced by a lattice constant should be finite in 
general. That a mobile hole carries a long-range dipolar 
spin distortion has been also studied in Refs. p!7lll9|. 

By contrast, in the present case, the hole is sitting at 
one of the two opposite poles of a dipole, not at its center. 
The holon itself has a finite hopping integral in Hh and 
can hop around. Nonetheless it is confined to an anti- 
meron by a logarithmic potential in Eq. H28|l . On the 
other hand, an anti-meron itself as a semi-classical vortex 
formed by the condensed spinous is immobile: it cannot 
propagate in terms of Hs- In other words, its effective 
mass is infinite, which is in contrast to a holon which 
has a finite hopping integral, or mass, according to Hh- 
Therefore, the hole dipole (both Type A and B) as a 
whole must remain localized or self-trapped in space. 

This distinction between the Shraiman-Sigga dipole 
and the present one can be traced back to a very singular 
effect in the t — J model, which has been inappropriately 
treated in the semi-classical approach of Shraiman-Sigga. 
In the following we further clarify this issue. 

According to Eq. (|22ll . there is a branch-cut connect- 
ing two poles (at respectively) of a hole dipole at 
the origin, as illustrated in Fig. 3. Such a branch-cut 
can be physically understood in terms of the irreparable 
phase string effect, which was discovered 01 as a generic 
property of the t — J model. 

A phase string is basically a sequence of -I- and — signs 
shown in [also see Fig. 1(a)], which is left by a spin- 
less hole on the path as it moves on a spin background 
described by the t — J model. The origin of the phase 
string is due to the disordered Marshall sign caused by 
the hole hopping Q . Physically it may be interpreted as 
describing displaced spins, by hopping, along the trans- 
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verse directions (in x-y plane). Specifically, the sign + 
and — in J^l distinguish the "backflow" spin f and J, ex- 
changed with a hole at each step of hopping. Since the su- 
perexchange term "respects" the Marshall sign, the phase 
string, describing the displaced Marshall sign, cannot be 
"repaired" through the superexchange term. Therefore, 
if one draws a circle around either the hole at B or the 
starting point at A, it will always intercept the phase 
string once (or odd times) [see Fig. 1(a)]. 

Pictorially, one may consider a bare hole created by 
the electron c-operator at a point A and then imagine 
that the spinless hole (holon) hops to the point B. It 
results in a dipole configuration with the two poles as its 
starting point A and ending point B, connected by the 
phase string in between [Fig. 1(a)]. Since the holon can 
reach B from different paths originated at A, the singular 
phase string is then replaced by (or relaxes to) a smoother 
dipole configuration. Nonetheless, the topology remains 
as represented by the branch-cut connecting two poles of 
the dipole as illustrated by Fig. 3. 

In a Neel state with the magnetization lying in x-y 
plane, the displacement of staggered magnetization mo- 
ments will reach the maximum, namely, tt phase per site 
(inversion of the spin polarization) on the hole path. It 
implies that going through the circle once, around the 
point A or B in Fig. 1(a), one finds the total spin dis- 
placement relative to the original Neel spin background 
is TT, modulo 2tt, which has the same topology as a meron 
or anti-meron discussed above. 

Finally, the localization of a hole dipole can be physi- 
cally understood in terms of the phase string picture too 
U^. Since the phase string ^ cannot be "self- healed" 
through, say, spin dynamics, a holon moves far away 
from its origin A will eventually lose its phase coher- 
ence as the product 1^ becomes uncertain due to the 
zero-point spin fluctuations in the Neel background. It 
means the propagating amplitude vanishes such that the 
holon must be localized or self-trapped around the origin 
A within a finite distance. The single electron spectral 
function for the one-hole doped case has been studied 
based on such a holon localization description, in which 
the dispersion of the low-lying peak is originated from 
the "spinon" sepectrum instead. The calculated spectral 
function well explains the corresponding photoemission 
experiments [Toj . 

Delocalization of holons will occur at higher doping 
when the spin background is changed (self-consistently) 
by the phase string effect. Such a holon deconfined phase 
is gene rally a superconducting state as discussed in Refs. 
fH Il2i| On the other hand, at low doping, without the 
destruction of spin ordering, delocalization of holons can 
only be realized in the reduced dimension when holons 
follow each other's paths such that the induced phase 
string frustrations can be eliminated. Such a metallic 
stripe in the spinon condensed phase will be discussed in 
the following section. 



IV. STRIPE INSTABILITY 

As discussed above, a hole-dipole object does not have 
a kinetic energy, or, its self-energy does not have a dis- 
persion, since it is self-trapped in space with an infi- 
nite mass. Thus the long-range dipole-dipole interaction 
among dipoles may become crucial in determining the 
low-energy physics in such a phase where kinetic energy 
is strongly frustrated. In this part we will discuss the 
origin of stripes based on this picture. 



A. Dipole-dipole interaction 

For a dipole located at the origin, Ea. H27l) can be re- 
duced to 



(zxr)(r-efe) 

2^ - Z X Gfc 



for |r| >> |efc| and a singular part 

V/ ~ -27r(z X ek)6{r) 



(35) 



(36) 



at |r| ^0 (|efc| ^0). 

Then, in terms of (|26|l . the dipole-dipole interaction 
between fc-th and fc'-th dipoles, located at and r^^/ , re- 
spectively, which are well-separated in space by jr^.^,' | = 
[r/c — Tfc' I 3> [e/cl, can be expressed in leading order ap- 
proximation by 



kk' 



d-d 



T s 2 
'JsPca TT 



2Tr Jsp^a 



(z X efc) • V/' (rfc) - (fc ^ k') 

2 (^fe ' ^kk')i^k' ' ^kk'^ 



''fefc' I 



efc • e^^- 



'^kk' I 







Here V^^T has a typical form of magnetic dipole-dipole 
interaction, with denoting the dipole moment. 

Consider orientations of dipole moments which min- 
imize the energy of the dipole-dipole interaction Eq. 
(|37|l at a fixed distance |rfcfc/|. Define cos (pk — ■ 
Tfefc'/ [rfefc'l and cosv^^' = ej,/ ■ Vkk' / \ek'\\rkk' \ ■ We 
have 



Trd — d 

^kk' 



iTTJsPlo? [efc I \ek' 



cos {(pik+^k-): 



|rfcfc' 

which is minimized at 



with 



-iT-d—d ^ 

min 



-^k' 



modulo 27r 



2-KJspla^ |efc| [efc/ 

Vkk'\^ 



(38) 



(39) 



(40) 



Other orientations of dipole moments always give rise 
to higher energies for a given \Tkk' \ ■ For instance, it is 
easy to see that for the case of random orientations < 
efc . efc, >= 0, {V^^") = 0. 
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B. Stripe instability 

Due to translational invariance, the center of a hole- 
dipole can be located at any spatial point with the 
same self-energy. Thus, to minimize H4()(l . the distance 
\^kk' I between two dipoles should be continuously shrunk. 
Without a balance from the kinetic energy, such collaps- 
ing seems an inevitable consequence of the dipole-dipole 
interaction. Of course, single hole-dipole picture is only 
meaningful at a scale jr^.^,' | > |efc| . Below such a scale, 
one has to take into account the detailed local physics to 
determine the hole-spin configuration and energy. 

It is important to note that this dipole-dipole "attrac- 
tion" does not simply result in a real-space pairing of 
holes. Due to the lack of kinetic energy cost, more dipoles 
would rather like to lump together, with orientations of 
dipole moments satisfying Eq. (|39l) , in order to lower the 
potential energy. They can thus either line up to form a 
stripe or phase separate into Wigner-crystal-like clusters. 
In reality, the existence of long-range Coulomb interac- 
tion may generally favor Q a stripe configuration over 
large Wigner-crystal-like clusters. Moreover, the latter 
possibility can be also considered as a further instabil- 
ity of stripes when they are attractive to each other and 
subsequently collapse into 2D clusters. Thus one can al- 
ways focus on single stripes first, based on the hole-dipole 
picture. In fact, we shall see below that, for instance, 
forming a metallic stripe will lead to a disassociation of 
individual dipoles and release additional kinetic energies 
which can further stabilize stripes against forming a 2D 
cluster. 

One possible loophole in the above argument is the ex- 
istence of some other forces which may be in competition 
with the dipole-dipole interaction. We first note that the 
dipole-dipole interaction could be quite weak. In a hole- 
dipole shown in Fig. 3, the holon sitting at one pole is 
actually not static: it generally will hop around the pole 
at A, in terms of Hh [Eq. 10)], such that {e^) may get 
strongly suppressed, and so does the dipole-dipole inter- 
action. It is a matter of competition between the hopping 
energy of holons (which contributes to the self-energy of a 
hole-dipole) and the dipole-dipole interaction in Eq. H38I) 
or Eq. H4UI) that self-consistently determines the actual 
dipole moment (e^) . 

In realistic situation, for example, a small amount of 
disorders can easily pin down hole-dipoles without the 
penalty from kinetic energy, such that the distribution 
of holes may still remain uniform on average. Based on 
a homogeneous picture, it has been previously shown 
that the AF long range state will persist over some finite 
doping concentration if interlayer spin-spin coupling is 
considered in the same model. The dipole-dipole inter- 
action will then be responsible for causing a cluster spin 
glass phase beyond some critical concentration of holes 
and eventually an insulating-superconducting transition 
takes place at a critical doping Xc ~ 0.043 

Even as a metastable state with local minima, such 
uniform phases may only evolve into stripe states of lower 



energy after a macroscopic time scale. So these uniform 
insulating phases are highly competitive with the inho- 
mogeneous phases, in a real experimental situation where 
a quenched disordered phase from high temperature may 
result in a spin glass phase 7, .20., 21, Tl^ .J3J , with the 
help of impurities, instead of a more ordered stripe phase 
at low temperature. 



C. Static stripes 

So the long-range dipole-dipole interaction will be re- 
sponsible for doped holes to "talk" to each other and for 
them to collectively collapse into stripes. It provides a 
mechanism for the instability of a homogeneity insulat- 
ing phase and the formation of stripes. But the local 
symmetry and hole-spin dynamics will ultimately decide 
which kind of stripes is most stable. In the following, 
we consider some general properties of stripes with hole- 
dipoles as elementary building blocks. 

Anti-phase. Let us consider a line-up of hole-dipoles 
along, say, i-axis with the dipole orientations satisfying 
the optimal condition Eq. Namely, {ek)^ = e^, 

(e^) = (— l)'^ey, where integer k is an index labelling 
dipoles. Here we assume that each dipole has the same 
size of a dipole moment |e| with a minimal spacing 2ex 
between two dipoles. For a site far away from the x-axis, 
the spin twist is then given by 



fe _ exVik - i-l)''eyXik 



E 



exVik 



k 



-sgn{yi), 



E 

k— — c<. 



{2exk + xo)^ + 



when \yi\ > |e| a. Here yik = yi and Xik = 2exk + xq 
according to the definition. 

Thus, an anti-phase (tt phase shift) is found across the 
stripe with 



Act) = 4>y>0 ~ <i>y<^ 



(41) 



Namely, a stripe composed of hole-dipoles is topologically 
always an anti-phase domain wall of spins. 

Metallic stripe. For Type A dipoles discussed in Sec. 
HI, a minimal choice of e is e^; = a and ey — 0, as shown 
in Fig. 4(a). If one simply line up these dipoles along 
the i-axis, then the corresponding stripe is an anti-phase 
domain wall with quarter-filled holons 



(42) 



(i.e., one hole per two sites), as illustrated in Fig. 5(a). 

Obviously the twists of spins due to the presence of 
hole-dipoles are all canceled out away from the stripe, 
such that the total potential energy is minimized. With- 
out spin frustrations, two stripes will not gain additional 
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Diagonal stripe from type B dipoles 

(b) 



FIG. 5: Two kinds of stripes as composed of two elementary 
types of liole-dipoles sliown in Fig. 4: (a) An anti-pliase, 
quarter-filling, metallic stripe; (b) An anti-phase, insulating 
stripe. 



energy by being close together and thus these stripes 
are also stable against a cluster formation. Furthermore, 
within each anti-phase stripe, the hole-dipoles will disas- 
sociate to gain additional kinetic energy. It is straightfor- 
ward to see now holons can move along the stripe with- 
out causing additional frustrations, just like in the one- 
dimensional (ID) t — J model jlj. Here the flowing of 
the "river of charges" along the stripe will gain energy 
~ — 0.9t/i for each holon at quarter-filling. This gaining 
of the kinetic energy is in contrast to the energy cost for 
a hole-dipole escapingaway from the stripe: f ~ Js 
according to Eq. (|28(l 2^ . The energy difference will sta- 
bilize such a metallic stripe of quarter-filling, at least in a 



local minimal sense. In literature, a metallic stripe based 
on the t — Jz model has been discussed in Ref. |2^] . The 
spectral density for a sin gle h ole in the stripe phase has 
been also studied in Ref. [2g. 

Based on the phase string picture, a metallic stripe 
discussed above can be regarded as a line which connects 
those branch-cuts between the starting and ending points 
of holons in Fig. 3. As emphasized before, a phase string 
is so singular that it can only be eliminated when one 
holon tightly tails another, which leads to a ID line of 
charges rather than a pairing of holes. Physically it is the 
phase string effect that ultimately threads holes together 
in a metallic stripe. 

As a natural extrapolation, one may picture a stripe 
liquid as composed of mesoscale stripes threaded by 
phase strings, after the melting of static stripes. Eventu- 
ally, say at high enough temperature, the most elemen- 
tary entities, i.e., hole-dipoles, should be also present, if 
spins are still ordered (Bose condensed) . Experimentally, 
how to detect such mesoscale stripes or even individual 
hole-dipole objects is an interesting subject to further 
explore in the present framework. 

Diagonal stripe. For a Type B dipole, its minimal 
dipole moment is shown in Fig. 4(b). Here one pole 
is located inside a plaquette while the pole of the holon 
is at a corner of the plaquette. Fig. 5(b) shows a nat- 
ural stripe configuration composed of the minimal Type 
B dipoles in Fig. 4(b). It is a stripe along the diago- 
nal direction, which remains as an anti-phase domain as 
the argument leading to Eq. H41|) still holds here. But 
it becomes an insulating stripe instead, as holons cannot 
directly move along the stripe. 

It has been argued that hole-dipoles of Type B may be 
more stable at low doping as it costs less superexchange 
energy. By forming a diagonal stripe, there is no obvious 
gain in the hole's kinetic energy, except for some virtual 
process. In a single hole-dipole case, such a virtual pro- 
cess may be understood as each holon continuously hops 
around and then back to its origin, leaving irreparable 
phase strings on the paths which result in a spin-dipole 
configuration. When these holes are lined up in the stripe 
configuration, the phase strings are canceled out each 
other during the virtual hopping process such that there 
is no spin twist left behind away from the stripe [Fig. 
5]. So for both types of stripe in Fig. 5, the healing of 
the spin AF order away from the domain wall indicates 
an effective kinetic energy gain for virtual hole hopping 
out of the domain wall. Furthermore, holons gains addi- 
tional kinetic energy in the metallic stripe, but the su- 
perexchange energy of the domain wall seems favored in 
the diagonal stripe, under a given number of holes. 

One may envisage some other kinds of static stripes. 
A detailed analysis of which type of stripes is most ener- 
getically favorable in a given doping concentration is be- 
yond the scope of the present paper, and will be explored 
elsewhere. Finally, we note that various charge order- 
ings due to hole-hole long-range (dipolar and Coulomb) 
interactions have been also studied numerically in Ref. 
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|27| . based on a spin-density- wave, or itinerant approach. 
Our strong-couling (t-J type model) approach seems to 
provide a consistent picture concerning the formation of 
static stripes, though many details still need a further 
comparison. 

V. DISCUSSIONS 

In this paper, we have studied how holes behave in a 
spin ordered background at low doping. In a doped AF 
Mott insulator described by a derivative model of the 
t — J Hamiltonian, known as the phase string model, the 
kinetic energy of doped holes gets strongly frustrated by 
the phase string effect. These holes induce nonlocal ef- 
fects in the spin background and behave like dipoles. In 
contrast to the hole dipole picture of Shraiman and Siggia 
[iGj based on the semiclassical treatment, a holon consti- 
tutes one of the two poles in the present configuration 
instead of sitting at the center of the dipole . Conse- 
quently, hole dipoles are localized or self-trapped objects 
in space rather than mobile ones in the Shraiman-Siggia's 
case. Such a localization or self-trapping of a hole can be 
intuitively understood directly in terms of the irreparable 
phase string effect [l9|. 

That the hole kinetic energy is strongly suppressed at 
low doping in the phase string model is in contrast to 
the slave-boson theories @ of the t — J model. In the 
latter cases, holes are relatively free moving on various 
RVB backgrounds. The difference arises from the fact 
that in the bosonic RVB background 8], short-range AF 
correlations are much stronger than in fermionic RVB 
theories, which are highly unfavorable to the hole motion. 

Due to the frustrations of the hole kinetic energy, 
the dipole-dipole interaction dominates the low-energy, 
long- wavelength physics. Holes correlate with each other 
through such an interaction which, instead of leading to 
a real-space pairing of holes, causes an inhomogeneous 
instability - stripe instability with hole-dipoles collaps- 
ing into a ID line. An anti-phase, quarter-filling metallic 
stripe can be regarded as a pile up of Type A dipoles, 
while an anti-phase, diagonal insulating stripe may be 
naturally considered as composed of Type B dipoles. 
Here the phase string effect provides the basic thread for 
forming stripes: holes tightly follow each other's paths 
in order to minimize the frustrations of the phase string 
effect and to release the kinetic energy. 

Theoretically, stripes and stripes related issues have 
been studied intensively in literature 0, motivated by 
experimental discovery of static stripes in curpates. 
In fact, insulating stripes as a mean- field solution was 
first obtained 29| [sO, |^ in the Hubbard and t — J 
models at low-doping even before the experiment. The 
present approach has gone beyond the mean-field approx- 
imation in dealing with the motion of doped holes in an 
antiferromagnet. It reveals that as a "renormalization" 
effect there does exist a residual long-range interaction 
for holes to "talk" to each other and to collectively form 



stripes. This approach provides a framework to further 
systematically investigate various properties associated 
with stripe phenomenon. Finally, in literature stripes as 
composed of ferromagnetic bonds beyond the t — J type 
model have been previously proposed js^ . But the origin 
of the underlying physics as well as the consequence are 
different. 

Some interesting questions remain to be answered. 
Firstly, if static stripes melt, are there dynamic stripes 
present at higher energy or higher temperature? If they 
exist, it is also a very challenge problem to mathemati- 
cally describe such dynamic stripes of "finite sizes" . Re- 
cently a Z2 gauge theory of "stripe fractionalization" has 
been proposed by Nussinov and Zaanen [Sl^ based on a 
"topological" order of the sublattice parity. The Z2 phase 
string depicted in Fig. 1(a) indeed originates from the 
change of the sublattice parity of spins due to the hole 
hopping. But the phase string model of Eqs. (O and 
is not simply a Z2 gauge theory because the original 
building blocks are holons and spinons, not in a stripe 
fractionalization description. 

Of course, in the end the holons are actually confined 
in the spinon condensed phase and elementary objects 
become hole-dipoles. A mesoscale stripe structure may 
be expected as remaining clusters at high energy or tem- 
perature, threaded by phase strings, which should even- 
tually break into hole-dipoles. But we note that the 
above dipole picture and stripe instability are all dis- 
cussed when spinons are Bose condensed. When such a 
condensation is gone at high temperature or larger dop- 
ing, the dipole and stripe pictures are no longer mean- 
ingful in this framework. 

Secondly, what is the relation between superconduc- 
tivity and stripe phases? So far we have not explored 
the possibility of superconducting condensation in the 
spinon condensed phase, where stripes present. In the 
phase string model, superconducting condensation hap- 
pens when holons are Bose condensed 0, It can 
occur at hig her doping where no spinon condensation is 
present [iJI- Thus, superconducting phase and stripe 
phase seem not necessarily concomitant in this frame- 
work. Nevertheless, the origin of superconductivity and 
that of stripes can all be attributed to kinetic- energy- 
driven mechanisms. Since the phase string effect plays 
the concrete role of kinetic-energy frustrations, various 
low-temperature phases including stripes, superconduc- 
tivity, and even the pseudogap phenomenon may be re- 
garded as different ways in reducing kinetic energies in 
such a system 0. 
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